Spectra of pinned charge density waves with background current 
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We develop techniques which allow us to calculate the spectra of pinned charge density waves 
with background current. 



I. INTRODUCTION 

The problem of a quantum mechanical particle whose 
motion is governed by a random time-independent 
Hamiltonian has been discussed in the literature for sev- 
eral decades. Because of its relevance to the electrons' 
motion in disordered conductors, the Hamiltonian which 
consists of a kinetic energy and a random scalar poten- 
tial is probably the most famous example of such a prob- 
lem. It is now well established that some of the wave 
functions of such Hamiltonians are localized, the phe- 
nomenon usually referred to as Anderson localization . 
A wealth of new types of behavior was discovered once 
random Hamiltonians constrained by certain symmetries 
were studied. For example, the most general random 
Hamiltonian would be complex hermitian. Constraining 
it to be real (in other words, imposing time reversal in- 
variance) changes the localization properties of its wave 
functions. An electron moving in a disordered conductor 
is described by a real random Hamiltonian. Turning on 
a magnetic field makes the Hamiltonian hermitian. This 
made the study of the crossover between real and hermi- 
tian random Hamiltonians easily accessible experimen- 
tally, and a subject of intense theoretical investigations 
0. 

Furthermore, it was discovered that a number of other 
symmetries can also be imposed on random Hamiltoni- 
ans, which change their behavior yet again. Altogether, 
there are ten symmetry classes of random Hamiltonians, 
distinguished by nine different constraints imposed on 
them [10. 

Yet other types of random Hamiltonians arise when 
they are constrained not by a symmetry but by certain 
requirements their spectra must satisfy. The most promi- 
nent example of that would be systems with bosonic exci- 
tations 0. To be specific, consider an energy functional 
E[<fi(x)], where <fr(x) is some function, which has the form 



E\ 



dx 



1 f d4{x) 

2 V dx 



h(<j),x) 



(1) 



Here h(<j)(x), x) is a random function whose form will be 
specified below. Suppose <f>o(x) is a minimum (local or 
global) of this functional. We would like to study the 
normal modes tf> n of oscillations around that minimum, 
described by the equation 



Eq. @ is a random Schrodinger-like equation. Yet it 
is not equivalent to a particle moving in an arbitrary 
random potential. d^h(<f>o(x), x) is a random function of 
x, but it is not an arbitrary random function. It is clear, 
for example, that e„ > for all n, which tells us that 
d^h((/)Q(x),x), although random, has to be constrained 
in such a way that all its eigenvalues are positive. 

In most applications of Eq. (JIJ, h(<p, x) is chosen to be 
a smooth function of <p and a rough function of x. For 
example, h((j),x) = A(x) cos(</> — x( x ))? where A(x) and 
x(x) are random functions of x uncorrelated at different 
x, and x( x ) i s uniformly distributed between and 2ir. 
This leads to the two point correlation function 

( h(<j>, x) h{<j>' , x')) = a cos(0 — </>') S(x — x'). 

Under this definition of h(cj>, x), the problem described 
by Eqs. and J3J) has been extensively studied in the 
context of pinned charge density waves |(| ■ Knowing the 
spectrum of Eq. J2J) allows, for example, to calculate the 
AC conductance of the charge density wave. It was first 
deduced in Refs. and that the density of states p(e) 
of Eq. (J2J) is given by p(e) = if e <C e c and if 4>o(x) is a 
global minimum of the energy functional Eq. Jy, where 
e c ~ cv 2 / 3 is the crossover scale. If, on the other hand, 
4>a(x) is but a local minimum of Eq. then p(e) ~ e, 
e < e c . 

A more detailed approach to the problem specified 
by Eq. (J2J was developed in Ref. 9]. It was shown 
in that work that the potential d^h(4>o(x),x) of the 
"Schrodinger" equation Eq. can always be represented 
as 

d 2 M<M^) = ^ + u 2 (*), 

where V(x) is some new random function. As a result, 
Eq. J2Jl is equivalent to 



where 



H 



n = 



1pn(x) 
(f>n(x) 



& + V(x) 

-£ + v(x) o 



(3) 



(4) 



dx 2 



ip n (x) + dlh((j)o(x),x) ip n (x) = e n ip n . (2) 



and a; 2 = e„. Now H. is an example of a Hamiltonian 
constrained by a symmetry and can be solved by tech- 
niques developed in that context. The symmetry of TL 
is usually referred to as chiral symmetry. It is expressed 
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by the relation a^Tia^ = —H, which holds true for any 
random V(x) Q. Here 173 is the usual Pauli matrix. 

The full solution of Eq. J3J), with the help of the map- 
ping to Eq. @, demonstrated that there indeed exists 
the crossover energy scale e c . At e n -€i e c all the eigen- 

functions ipn(x) are localized with the localization length 

_i 

£n ~ £c 2 , which is independent of n. At e„ ^> e c , the 
eigenfunctions are still localized, but with the localization 
length which increases with e„ as ~ e„. The density 
of states at e„ 3> e c is given by p(e) ~ . Finally, the 
density of states at e„ <§; e c is given by p(e) ~ e 5 if <j> (x) 
is a global minimum of Eq. (JIJ and by p(e) ~ e if </>o(a;) 
is a local minimum of Eq. Q . 

In this paper we would like to show that the same 
techniques which proved useful in solving Eq. J5J can 
also be used to solve another related problem, which we 
formulate below. Consider an equation 



J 



dx 2 



1< ~<L + d * h ^' x ^ 



0. 



(5) 



where j and 7 are some parameters. If 7 = 0, then 
this equation is equivalent to the minimization condition 
of the energy Eq. If, on the other hand, 7 > 0, 

then this defines a new problem. Eq. (JSJ) was suggested 
in Ref. .10] to describe pinned charge density waves with 
background current. Normal modes of oscillations of such 
pinned charge density wave are given by 



-3- 



~ f — + dlh{cj> (s {x),x) 



dx 2 ' dx 

(6) 

where <fio (x) is a solution to Eq. (JSJ) ■ In this paper we are 
going to present the solution to the problem defined by 
Eq. © and Eq. ®. 

The operator in the square brackets of Eq. 10 is non- 
hermitian. As a result, the eigenvalues e„ do not have to 
be real. It is well known, however, that the eigenvalues of 
equations of the type Eq. © always lie along one dimen- 
sional curves in the complex plane In this paper, we 
show that these curves take the shape depicted on Fig. 
^ The fork point tf cannot be found exactly. However, 
it is still possible to define the crossover scale e c . All the 
states with energy less than e c have the same localization 
length £. As the energy is increased past e c , the localiza- 
tion starts to grow and eventually diverges at ef (which 
is always bigger than e c ). The states corresponding to 
complex values of energy are delocalized. 

We introduce the technique which allows to compute 
e c and £. We calculate them in the regime where 
7/(ja) 1 / 3 3> 1 and find them to be e c = 7 2 /4j and 
£ = 2j 2 /a. We also calculate the density of states along 
the part of the curve depicted on Fig. ^ which lies on the 
real axis and find it to be p(e) ~ e. 

The rest of this paper is organized as follows. In section 
|h|we map the problem defined by Eq. © to a random 
chiral Hamiltonian. In section 1TTT1 the localization length 
of normal mode oscillations will be derived in the limit of 
large background current and in section llVl this result will 




FIG. 1: The eigenvalues e n of the Eq. lj|]J lie on this curve in 
the complex plane. 



be compared with the usual Larkin length of the prob- 
lem. Section IV1 contains the calculation of the density of 
the low lying states. Finally, in Section I VII we discuss 
the response of a pinned charge density wave with back- 
ground current to an electric field using the formalism 
developed here. 



II. MAPPING TO A CHIRAL HAMILTONIAN 

The first step of the solution is to map Eq. © to a more 
tractable form of a random chiral Hamiltonian. Follow- 
ing the general guidelines of Ref. |9j] , we treat Eq. © as 
the equation of motion of a particle with coordinate </>, 
moving in time x, in the presence of viscous forces given 
by the 7 term. We then pass from the Lagrangian to Eu- 
lerian description of the motion. Introducing the veloc- 
ity function u(4>o(x)) = d x 4>o{x) the equation of motion 
Eq. © becomes 



jd x u + jud^u + ju = d^,h(4>o(x),x). 



(7) 



The equation thus obtained is similar to the equation 
of motion of a one dimensional fluid without pressure, 
usually referred to as Burgers equation. However, it also 
features the 7U term, which can be interpreted as a kind 
of viscous friction. The Burgers equation with this term 
is novel and has not been considered in the literature 
before. 

Performing a derivative on Eq. Q with respect to <f> 
results in 

jd x d^u + j(d<f,u) 2 + judlu + ^fd^u = d\h. (8) 

We define the gradient of the velocity along the solution 
to Eq. ©, 0o, as F(x) = d^u((j)o{x),x). F(x) can be 
related to h((j), x) with the help of Eq. (|SJ 



dF 



j-j 1 - ]F + jF = d^h{4> Q {x),x). 



dx 

Thus Eq. jnj becomes 



(9) 



d 2 id dF . 7 _' 

1 \-F + —F 

dx 2 j dx dx j 



J 
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or in a more symmetric form 



-r + F + ^r 

dx j _ 



d ' 

dx 



?Pn = —Ipn- 
J 



(11) 



The operator on the left-hand-side of Eq. is not Hcr- 
mitian. This may be remedied by using a trick due to 
Hatano and Nelson 01 ■ Writing 



if>n(x) 



the equation becomes 



d p 7 
dx 2j 



d p 7 
dx 2j 



(12) 



-^ n .(13) 
J 



Thus we mapped our equation into an eigenvalue prob- 
lem for a Hermitian operator. This operator can also be 
rewritten in the form similar to Eq. Ip}. 



H 



4>n(x) 



I0n_ ( 1pn(x) 
VJ \4>n{x) 



where Ti is still given by Eq. uj^ = e„, and 



V{x) = F{x) + 



2j' 



However, as argued by Hatano and Nelson, the transfor- 
mation Eq. (fT^|) is valid only as long as ip n decays asymp- 
totically as exp(— with the localization length 
£ < 2,7/7. F° r V'n whose localization length obeys this 
condition, Eq. l]l l|l is equivalent to Eq. (|13J) . and there- 
fore, e n is real and positive. For other eigenfunctions ip n 
whose localization length is larger than 2j /j, Eq. I|13|) is 
no longer equivalent to Eq. i|ll[l . 

Notice that the localization length £„ of ip n is related 
to the localization length £ n of ip n as in 



6. 



as long as £ < 2,7/7. 



2_ 

2j 



(14) 



III. THE LOCALIZATION LENGTH 

Random chiral Schrodinger equations of the form i|13|) 
have been investigated by A. Comtet, J. Debois, and C. 
Monthus 0]. It was determined that the important pa- 
rameter for these equations is w c = ^fj ( V[x) ), and cor- 
respondingly e c = Wj. For all wavefunctions whose en- 
ergy io n < lu c , the localization length is constant and 
is given by £ = 1/ (V(x)). For w„ > oj c , the localiza- 
tion length quickly increases, being asymptotically pro- 
portional to u>\. 

The crucial test for our theory is, therefore, whether 
( V(x) ) is bigger or smaller than 7/2,7'. As we will see be- 
low, (V(x)) > -f/2j. As a result, Eq. (fT*5|) is equivalent 



to Eq. (|llf> and consequently to Eq. ©, at e < e c . Thus 
at e < e c all the eigenvalues of Eq. © are real and posi- 
tive. On the other hand, there exist some e/ > e c where 
the localization length of ipn becomes equal to 2j/-f. At 
that point, the wave functions ipn become delocalized, in 
accordance with Eq. (|14|l . At Re e > e/, the eigenval- 
ues of Eq. © are no longer real and come in complex 
conjugate pairs. This justifies the picture presented on 

Fig.m 

In what follows, we proceed to calculate e c . In terms 
of V(x), Eq. © becomes 



dx 4j 2 



dlh{<j) (x),x) 



(15) 



This equation has the form of a Langevin equation with 
random force given by d^h, whose correlator is given by 
(see Ref. Q for justification of dih as a random white 



(dlh((j} (x),x) 9|/i(>o(y),j/)) = a S(x-y). 
In accordance with the theory of Langevin equations, 

if 

f + 9 (V) = M 

with f(x) being uncorrelated at different values of x, 
(f(x)f(x')} = a5(x — x'), then the probability P(v, x) = 
(6(v — V(x)) of observing V(x) = v at the position x 
obeys the Fokker-Planck equation 



dP 
~dt 



d_ 

dv 



ad . . 
2d~ V +giv) 



P(v,x). 



In the present case, we are interested in the probability 
of observing a particular value for V(x) together with the 
probability that a particular solution (f)o(x) of Eq. © 
is chosen. The relevant quantity describing this joint 
probability is 



V{v,x) 



5{v-V{x)) 
p(x) 



where p(x) is the density of solutions (po{x) which in turn 
obeys the continuity equation 



dp 

dx 



+ pV = 0. 



As a result, we find the Fokker-Planck equation 

,,2 



dP _ d_ 

dx dv 



a d 
2p dv 



4j 2 



V + \vV. 



Here A = 1, however, we will keep the more general no- 
tation of A for convenience at a later point in the calcu- 
lations. 
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A common approach to the Fokker-Planch equations 
is to map them into the Schrodinger equation with the 
help of 



This gives 

Here 

U(v) = 



V = exp 



dx 



v 3 j 2 
3a 



2j2 Qy2 



vy 
4a 



U(v) 



7 2 T 2 

J 4 / 2 / -i i \ \ 
2S° "to" ~ W(1 + A) 



(16) 



(17) 



The Feynman path integral formulation of the 
Schrodinger equation Eq. (|16|) allows us to find the aver- 
age of V using 



£(V) = 



JW(x) e- s Wf*dxV(x) 
JVV(x)e- s l y ] 

4r log / T>V(x)e- 8 M 
d\ J 



(18) 



for £ — > oo. Here S^V] is the imaginary time action which 
corresponds to the quantum mechanics Eq. . 



5[V] 



dx 



J 



dV 



2a V dx 



+ U(V) 



(19) 



At large I the path integral Eq. I|18fl is dominated by 
its ground state energy. The average V may now be 
calculated as 



(V) 



dX 



(20) 



A=l 



where Eq is the quantum mechanical ground state energy 
for a particle whose classical action is given by Eq. 119(1 
and whose Schrodinger equation reads 



'2fd^ 



U(v) 



If 7 = 0, then the ground state energy, together 
with its derivative with respect to A, can be estimated 
simply from dimensional analysis as a 1 ^ 3 /j 2 ^ 3 , to give 
uj c ~ a 1 ' 3 /] 1 ^. This agrees with Ref. 9). If 7 > 0, 
then dimensional analysis is of no help, since lo c can now 
depend on the dimensionless ratio 7/C/a) 1 / 3 . The only 
way to find oj c and thus the localization length of the low 
lying states is by finding the ground state energy Eq. 

By a suitable rescaling v = yo 1 / 3 / j 2 / 3 , E = 
we can bring the Schrodinger equation to 

the form 



^ = EV, (21) 



where 



2(a 3 Y 



In general, it is not possible to find the ground state 
energy of this Schrodinger equation exactly. We are going 
to find it only in the limit when c ^> 1. The potential in 
Eq. I|21|) has two minima. The minimum at positive x is 
the global minimum and is located at 



(22) 



The simplest approximation to the ground state energy 
is to set it equal to the value of the potential at the 
minimum. This is given by 

£0,0 = -(1 + X)c - + 0(c- 5 ). 

To add quantum fluctuations to the problem, we have to 
consider the Schrodinger equation (|21|l evaluated around 
the minimum Eq. I|22l) . Approximating the potential of 
Eq. (|21|1 by a quadratic potential around the point x m - ln 
we find the oscillator ground state energy to be 

^/= C+ ^^+0( C - 5 ). 

In principle, the cubic and quartic term in the expansion 
of the potential around the minimum also contribute to 
the ground state energy, but these only have dependence 
upon A at the 0(c~ 5 ) level and thus do not contribute to 
(V) at lower orders. 

Gathering the terms, it is seen that as a function of A 



(V) 



dEo 
' dX 



Eo,o + Eo. q f 



dX 



J- 



Finally, substituting A = 1, we find 

^ = 27 + ^ + ° ( ^ 5) - 

Therefore, at large 7 the energy scale lu c is given by 

7 



(23) 



We see that (V) > j/2j, and the localization length of 
the low lying hermitian states ip n is consequently smaller 
than 2^/7. This allows us to deduce that Eq. © is indeed 
equivalent to Eq. (|13fl at e n < e c . The localization length 
of ip n at e„ < e c can be deduced from Eq. 114(1 to be 



Z=[(V)- 



2 7 2 



+ 0( 1 - 1 ). (24) 



The calculation presented here justifies Fig. ^ at large 
7. In order to see that the picture of low lying localized 
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FIG. 2: { y ) in the ground state of Eq. (12 1 1 as a function of 
c. The dots represent { y ) evaluated numerically at various 
values of c. The solid line is ( y ) calculated in the harmonic 
approximation to the potential in Eq. (1211 . The dash-dotted 
line represents the first two terms of the perturbative calcu- 
lation c + . Finally, the dashed line is a straight line of the 
slope 1, which demonstrates that (y) > c. 



states persist at all values of 7 we need to show that ( V ) 
is always greater than "f/2j. This amounts to showing 
that in the Schrodinger equation Eq. (|21[l the following 
relation always holds true 



, y ; 



dy \My)\ 2 y = - ^ 



> c. 



(25) 



A=l 



Here ipo (%) is the ground state wave function of Eq. (|21[1 , 
at A = 1 . We do not know how to show this analytically. 
We investigated Eq. I|21|) numerically. Fig. [21 shows the 
values of — dE /d\ = ( y ), at A = 1, plotted versus c. We 
see that Eq. I|25|l does seem to hold. 



IV. THE LARKIN LENGTH 

In the analysis of Eq. JSJ it is customary to introduce 
the notion of Larkin length L. L is defined as size of 
a box in which the average (<p 2 ) becomes of the order 
of 1. Thus at distances bigger than L, 4>o{x) becomes 
rough, while at shorter distances it can be thought of as 
smooth. We are now going to see that in our problem 
the localization length of low lying states of Eq. © is of 
the order of Larkin length. 

To calculate the Larkin length, we follow the procedure 
described in Ref. 01 • We consider the equation Eq. JSJ. 
Let us assume that </> is small and neglect the </> depen- 
dence of h. In Fourier space this equation of motion is 



Thus 



Ifo 



j 2 k 2 +j 2 k 4 



Now the Larkin length is defined as the value of L for 
which 



1 



dk 



a 



- fk 4 



, -1 2tt7 2 /c 2 

La ja i . . . 
• tan-^/j). 



V 



From the asymptotic form of tan 1 in the large-7 limit 
it is found that 



L 



21 

a 



:31 



•0( 7 - 



At leading order L turns out to be proportional to the 
localization length obtained above Eq. H24|) . This co- 
incides with the behavior of the pinned charge density 
waves without the background current [9(. However, the 
localization length and the Larkin length have different 
functional dependence on 7 in lower orders. 



V. DENSITY OF STATES 

In this section we calculate the density of low lying 
states for our problem. Due to the form of the potential 
Eq. l(T7|) which does not contain a cubic term, the density 
of low lying states may be easily calculated using the 
methods of ref. |5J|. This is not immediately obvious, so 
in this section we will repeat the argument. 

At energies below an energy of the order of (V) the 
energy eigenvalues are real. Above, the eigenvalues trace 
out a one-dimensional spectrum in the complex plane, 
Fig. ^ Let us concentrate on low lying states e„ < e c , 
which are all real. To calculate their density of states 
p(e) we use the equivalence of Eq. (jBJ) and Eq. I|13|l in 
this regime. 

The integrated density of low- lying states can be found 
using the discussion in Ref. [5j as 



N{u)) = {6 yj dx V(x)+a 
a = -log(w), 



with the average being over realizations of V(x). The 
density of states is then obtained from p(uj) = dN/dui. 
Through a standard representation of the (5-function 



N{u) 



da ia( /*2 dxV(x)+aj 



2tt 



and defining S(a) by 

« P ,-«,o)) = ( ffi p( fa £W(*))), <*) 
with £ = X2 — x\ the integrated density of states is 
N(oj) = r ^ e ««-«(a)_ 

^ ' J '-co 2^ 

At large I , S(a) is independent of I as will be argued 
below. We assume in the following that a certain value 
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of £ maximizes the probability of observing a fluctuation 
of size — logo;. This is justified, since due to the form 
of the potential U(V), Eq. ((17(1 . it is not very likely to 
observe a very large negative fluctuation, and it is more 
likely to see a smaller negative fluctuation, which exists 
over a longer interval of x. We now approximate N{ui) 
by its value at the saddle point, determined from 

ia-l^- = 0, 
oa 

where a = ao(£) at the saddle point. Thus N(u>) ~ 
exp(iaao— IS (ao(£))) and the maximization with respect 
to £ results in 



= ia ^-S(a (£)) 
= -S(a (£)). 



£ 



dS_da 
~dcVd£ 



Defining (3q = — iaQ it is seen that N(u>) ~ luq where (3q 
should be chosen such that S(f3 ) vanishes. 

Making use of our results from section IIIII where the 
Schrodinger equation Eq. ((16(1 was found with the poten- 
tial taking the form Eq. (|17f) . Eq. I(26[) becomes 



J VV{x) exp (S[V}) exp (-#, Q datV) 
fT>V(x) exp(-S[V\) 



,(27) 



with S[V) given by Eq. (JTJJJ. At large lengths £ the path 
integrals are dominated by the ground state energy, i.e. 



exp(-£S{fto)) 



exp(-tE {-l-A + /? }) 
exp(—£E {—1 — A}) : 



where the notation Eq {s} represents the ground state 
of the Schrodinger equation Eq. 1)16(1 with the potential 
Eq. I|17fl where s is substituted in place of 1 + A. Fur- 
thermore, since all terms except the linear are invariant 
under V(x) — » —V(x), we see that Eq { — 1 — A + /3q} = 
E {1 + A — /? }. Thus (3 = 2 + A obviously solves the re- 
quirement. Note that the above argument hinged on the 
fact that the potential U(V) does not contain a cubic 
term. We conclude that with A = 1 



and 



p(u) 



dN(e) 
de 



(28) 



We would like to compare this with the result obtained 
for the pinned charge density wave problem with 7 = 0. 



There p(u>) 



if 0o is a global minimum of the en- 
; 3 if the minimum 



ergy functional Eq. and p(oo) 
is local. At 7 > 0, however, Eq. (jSJ) is no longer a min- 
imization condition of any functional and the notion of 
global minimum no longer exist. 



VI. DRIVEN PINNED CHARGE DENSITY 
WAVES WITH BACKGROUND CURRENT 

Consider a pinned charge density wave driven by an 
electric field at frequency uiq (see e.g. Ref. [l4j for a 
review) , 



d<t> 
'~dt 



J 



dx 2 



1-T- + d^h(<f>, x) = E(x) cos(woi). 



To find the linear response to a small electric field, we 
write 4> — 00 (%) +'4 > (x,t), where 4>o{x) satisfies Eq. JSJ, 
and find 



d 



ip(x,t) = E(x)e luJot . 



The solution to this equation can now easily be found 



(x,t) = 4> Q {x) + e 



dy E(y) 



^ n {x)ij} n (y) 



ivujq + e r , 

IL 

where ip n and e„ are defined in Eq. ©. The current 
carried b; 
given by 



carried by the charge density wave is I ~ ^ , and is thus 



I{x) ~ iuiQe 



dy E(y) 



ijj n (x)ip n (y) 
iuu) + e„ 



The properties of t/j n , e n found in this paper help to deter- 
mine the response of the charge density wave to a small 
external electric field. For example, if the electric field 
E(x) is uniform in space, we can use the results of this 
paper which show that at e„ < e c , the wave functions 
ip n {x) are all localized with the same localization length 
£. Therefore, J dx ip n (x) oc and is independent of n, 
as long as e„ < e c . On the other hand, the wave func- 
tions which correspond to large Re e„, are delocalized 
and oscillate fast, so that J dx ip n {x) goes to zero with 
increasing n quickly. Therefore, the current reduces to 



iu> t 



de p(e)- 



1 



iui Q t 



ivujQ log 



IVUJQ 



e 



IVUlQ 



(29) 



where p(e) ~ e, as found in this paper, Eq. 1)28(1 . 



VII. CONCLUSIONS 

We have demonstrated how the problem of normal 
mode oscillations of the one dimensional pinned charge 
density waves with background current may be solved 
by mapping into chiral random Hamiltonians. Using the 
methods of Comtet et al 01 we have determined the lo- 
calization length £ of the low lying normal modes. This 
localization length turns out to be proportional to the 
Larkin length of this system. The density of the lowest 



7 



lying states has also been obtained, giving the power law 
p(ui) ~ w 3 . This result differs from the result p(u>) ~ uj a 
of Refs. 0,0] due to the fact that turning on background 
current removes the notion of a global ground state. 
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